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ABSTRACT 

We developed a non-linear Monte Carlo code based on the LP method proposed by Stern et 
al (1995). We use it to simulate an extended non-thermal X/^y ray source coupled with a cold 
accretion disk. The source radiates by inverse Compton process. We show that very low optical 
depths (T ~ 5.10 - ) are required for the 2nd Compton order to be negligible. Such optical 
depths can be reached only if the source is very close to the disk (corona-like geometries). And 
even in this configuration, 2nd Compton order is negligible only for the hardest observed spectral 
slopes. We also account for Compton reflection on the disk. This kind of coupling leads to a 
strongly angle dependent spectrum with a reflection component much stronger than in standard 
isotropic models, leading to unobserved reflection dominated spectra at low inclination. 
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1. INTRODUCTION 

The high energy spectra of Seyfert galaxies and black hole binaries in their soft 
hard state are very similar. Their basic characteristics are an underlying powerlaw 
with photon index a ~ 1.6 — 2., a Compton reflection component with an Fe Ka 
line and a high energy cutoff above ~ 200 keV. The powerlaw is generally though 
to originate through Comptonisation of soft photons by a population of thermal 
hot electrons. Here we will focus on non-thermal models following Henri & Petrucci 
(1997, HP97) who proposed a model where the hard X-rays are emitted by a non- 
thermal optically thin point source. This model reproduces the basic features of 
the high energy spectra, but requires high inclination angle for all systems unless 
one assumes that a fair part of the disk is highly ionised (Malzac et al. 1998). 
However considering a point source, as done in those works, is unrealistic (HP97) 
As the source was assumed to be optically thin the second Compton scattering order 
(CSO) was neglected and the calculations performed in the Thomson limit. Here 
we developed a Non Linear Monte Carlo Code (inspired from Stern et al., 1995) to 
investigate the effects of source extension taking properly into account the radiative 
transfer and the Klein Nishina effects. 
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2. MODEL ASSUMPTIONS 

The X-ray source is modeled by a cylinder at some height above a disk. This ge- 
ometry could be physically realized by a strong shock terminating an aborted jet 
(Henri & Petrucci 1997). The source radiates by inverse Compton process illu- 
minating the disk. The seed photons are provided by the thermal disk emission 
at fixed blackbody temperature kT r . This disk is represented by an infinite slab 
which radiates only the reprocessed energy from the hot source without internal 
dissipation. In addition a fraction of the high energy photons intercepting the disc 
is Compton reflected before being absorbed leading to the formation of the so called 
reflection component. The relevant geometric parameters characterizing the source 
are x an d (see Fig. 1). The cylinder contains relativistic electrons or pairs with a 
fixed energy distribution: 



These leptons cool on the soft photon field. We assume that some acceleration 
process is active in the X-ray source to keep the above distribution fixed. The 
equilibrium cylinder optical depth depends only on the geometry and electron dis- 
tribution and remains to be determined self-consistently. The results are scale and 
luminosity independent since we neglected photon-photon pair production and pair 
annihilation. 

The distribution parameters are constrained by the observations: 2 < s < 3, 
7o < 200, kT r > leV. j m in < 10 otherwise the lack of low energy electrons would 
cause an unobserved break above 1 keV in the photon spectra. The most drastic 
contraint arises from the requirement of the second and higher CSOs to be negligible. 
Indeed, due to the high mean electron energy, the combination of CSOs leads to a 
photon spectrum extending well above lMeV. Thus the model requires very small 
optical depths typically r ~ 5.10 -4 (green and yellow region of Fig. 2). The relative 
strengh of the 2nd CSO depends also on s, 7 m i„ and the viewing angle (Fig. 5,6). 
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3. RESULTS 

Our NLMC code computes both the optical depth and the resulting photon 
spectra. This code has been tested successfully against a linear code (Podzniakov 
et al, 1983). In this version of the code, the power supplied to the electrons as well 
as their distribution are fixed. At each time step the power is supplied in the form of 
an injection of electrons. In balance, the number of electrons is decreased according 
to the Compton losses. So that, the optical depth 'evolves' until it reaches an 
equilibrium. We performed over 100 simulations covering a wide range of parameter 
values. The optical depth can be defined along the height of the cylinder or along 
its radius: 



where K is a function of 9, x, 7mm, 70j s > P and kTr. K is the lowest optical depth 
that can be achieved in any direction. 

Using our NLMC simulations we studied the dependence of K on these param- 
eters. While x has a slight influence on K, increasing K of at most a factor of 2 
between x — 10~ 3 and x = 10 3 , K is very sensible to the solid angle sustended by 
the cylinder as seen from the source. Fig. 3 shows K as a function of cosO K for 
X = 1) 7o = 50, J m in = 2, j3 = 2, s = 3, ktr — 50eV. K decreases strongly with 9, 
minimum values are obtained for 9 ~ n/2. 

kTr as pratically no influence on K provided it is lower than 0.1 keV. K is very 
sensible to the mean lorentz factor of the electron distribution. Fig. 4 shows K 
versus < 7 > for x — L cos ® = -83, kTr = bOeV and various values of 7 m m, 7o> P 
and s. Clearly, K is minimum for a slab coronna geometry (x = and cos9 = 0) 
and large values of < 7 > . We will now focus on an extreme case of corona geometry 
to study a configuration which provides the lowest achievable K values. We fixed 
70 = 200 (its maximal observational limit compatible with kTr = leV). Fig. 2 
shows the evolution of K as a function of "f m in for s = 2.1 (upper curve) and s = 3 
(lower curve) (x = 2. 10~ 3 , cos9 = 0, 70 = 200, ,3 = 3.). Note that in this limit of 
low Xi K ~ t h- In the lower region the 2nd CSO is negligible for s < 3.; in the 



m =-^a T H = K(l+ X ) 



tr = -^<t t R = K(1 + l/x) (2) 
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middle region it is negligible at least for s = 2. In the upper and darker region 
the 2nd CSO is never negligible for 2 < s. These regions have been determined 
qualitatively by simulations of a slab geometry with the linear code. The green 
region can not be reached even for the extreme value of "f m in = 10. However, the 
curve s = 2.1 crosses the middle region. We thus see that the model can only 
reproduce the spectra with the hardest observed slopes. 

Fig 5 shows the spectra obtained with s = 2.1, j m i n — 10 (th ~ 4 1CP 4 ) for 
different inclination angles. The 2nd CSO is well below the first order for E < 
1 MeV. It becomes dominant above 1 MeV, but at these energies the flux can 
be considerably lowered by photon-photon absorption. Fig. 6 shows the spectra 
obtained with s = 3, jmin = 10 (th ~ 8. 10~ 4 ) for different inclination angles. The 
2nd CSO is important for all viewing angles. The spectrum resulting in the sum of 
the two CSOs does not show any cut-off around 200 keV. These spectra show that 
the previous results concerning the angular dependence of the spectra (HP97) are 
not sensibly affected by the source extension. The reflection component appears a 
bit lower, relaxing somewhat the constraints on the inclination angle (Malzac et al., 
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FIGURE 6. s = 3 

1998). However the spectra are still reflection dominated at small inclinations. 
4. CONCLUSION 

Our self consistent simulations showed that low optical depths could be reached 
only if the source is very close to the cold disc (corona- like geometry). Moreover, 
even in this case, the 2nd CSO is negligible only for the hardest observed spectral 
slopes a ~ 1.6. Lower optical depths could be obtained if consequent dissipation 
occured in the disc. However, the source extension does not suppress the unobserved 
dominent reflection component at small inclinations. 
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